Biological efficacy of simulated radiolabeled Lipiodol® ultra-fluid and microspheres for various beta emitters: study based on VX2 tumors.

BACKGROUND: Radioembolization is one therapeutic option for the treatment of locally early-stage hepatocellular carcinoma. The aim of this study was to evaluate the distribution of Lipiodol® ultra-fluid and microspheres and to simulate their effectiveness with different beta emitters (90Y, 188Re, 32P, 166Ho, 131I, and 177Lu) on VX2 tumors implanted in the liver of 30 New Zealand rabbits. RESULTS: Twenty-three out of 30 rabbits had exploitable data: 14 in the group that received Lipiodol® ultra-fluid (group L), 6 in the group that received microspheres (group M), and 3 in the control group (group C). The histologic analysis showed that the Lipiodol® ultra-fluid distributes homogeneously in the tumor up to 12 days after injection. The X-ray µCT images showed that Lipiodol® ultra-fluid has a more distal penetration in the tumor than microspheres. The entropy (disorder of the system) in the L group was significantly higher than in the M group (4.06 vs 2.67, p = 0.01). Equivalent uniform biological effective doses (EUBED) for a tumor-absorbed dose of 100 Gy were greater in the L group but without statistical significance except for 177Lu (p = 0.03). The radionuclides ranking by EUBED (from high to low) was 90Y, 188Re, 32P, 166Ho, 131I, and 177Lu. CONCLUSIONS: This study showed a higher ability of Lipiodol® ultra-fluid to penetrate the tumor that translated into a higher EUBED. This study confirms 90Y as a good candidate for radioembolization, although 32P, 166Ho, and 188Re can achieve similar results.

Performance study of a 360° CZT camera for monitoring 177Lu-PSMA treatment.

BACKGROUND: The aim of this study was to investigate the quantification performance of a 360° CZT camera for 177Lu-based treatment monitoring. METHODS: Three phantoms with known 177Lu activity concentrations were acquired: (1) a uniform cylindrical phantom for calibration, (2) a NEMA IEC body phantom for analysis of different-sized spheres to optimise quantification parameters and (3) a phantom containing two large vials simulating organs at risk for tests. Four sets of reconstruction parameters were tested: (1) Scatter, (2) Scatter and Point Spread Function Recovery (PSFR), (3) PSFR only and (4) Penalised likelihood option and Scatter, varying the number of updates (iterations × subsets) with CT-based attenuation correction only. For each, activity concentration (ARC) and contrast recovery coefficients (CRC) were estimated as well as root mean square. Visualisation and quantification parameters were applied to reconstructed patient image data. RESULTS: Optimised quantification parameters were determined to be: CT-based attenuation correction, scatter correction, 12 iterations, 8 subsets and no filter. ARC, CRC and RMS results were dependant on the methodology used for calculations. Two different reconstruction parameters were recommended for visualisation and for quantification. 3D whole-body SPECT images were acquired and reconstructed for 177Lu-PSMA patients in 2-3 times faster than the time taken for a conventional gamma camera. CONCLUSION: Quantification of whole-body 3D images of patients treated with 177Lu-PSMA is feasible and an optimised set of parameters has been determined. This camera greatly reduces procedure time for whole-body SPECT.

International recommendations for personalised selective internal radiation therapy of primary and metastatic liver diseases with yttrium-90 resin microspheres.

PURPOSE: A multidisciplinary expert panel convened to formulate state-of-the-art recommendations for optimisation of selective internal radiation therapy (SIRT) with yttrium-90 (90Y)-resin microspheres. METHODS: A steering committee of 23 international experts representing all participating specialties formulated recommendations for SIRT with 90Y-resin microspheres activity prescription and post-treatment dosimetry, based on literature searches and the responses to a 61-question survey that was completed by 43 leading experts (including the steering committee members). The survey was validated by the steering committee and completed anonymously. In a face-to-face meeting, the results of the survey were presented and discussed. Recommendations were derived and level of agreement defined (strong agreement ≥ 80%, moderate agreement 50%-79%, no agreement ≤ 49%). RESULTS: Forty-seven recommendations were established, including guidance such as a multidisciplinary team should define treatment strategy and therapeutic intent (strong agreement); 3D imaging with CT and an angiography with cone-beam-CT, if available, and 99mTc-MAA SPECT/CT are recommended for extrahepatic/intrahepatic deposition assessment, treatment field definition and calculation of the 90Y-resin microspheres activity needed (moderate/strong agreement). A personalised approach, using dosimetry (partition model and/or voxel-based) is recommended for activity prescription, when either whole liver or selective, non-ablative or ablative SIRT is planned (strong agreement). A mean absorbed dose to non-tumoural liver of 40 Gy or less is considered safe (strong agreement). A minimum mean target-absorbed dose to tumour of 100-120 Gy is recommended for hepatocellular carcinoma, liver metastatic colorectal cancer and cholangiocarcinoma (moderate/strong agreement). Post-SIRT imaging for treatment verification with 90Y-PET/CT is recommended (strong agreement). Post-SIRT dosimetry is also recommended (strong agreement). CONCLUSION: Practitioners are encouraged to work towards adoption of these recommendations.

Prognostic impact of bone metastases detected by 18F-DOPA PET in patients with metastatic midgut neuroendocrine tumors.

OBJECTIVES: Bone metastases (BM) may influence negatively the prognosis of midgut neuroendocrine tumors (NET). The diagnostic sensitivity of 18F-DOPA PET for midgut NET and associated metastases is high. This study aimed to assess the prognostic impact of BM detected by 18F-DOPA PET in metastatic midgut NET. METHODS: All patients with a metastatic midgut NET, who underwent a 18F-DOPA PET between June 2011 and June 2018, were included. BM were defined following imaging criteria and were classified as poly-BM or oligo-BM, according to their number (< 5 or ≥ 5, respectively). The variables associated with the presence of BM were evaluated by logistic regression. The factors associated with overall survival were explored by Cox regression models. RESULTS: Among 155 patients included, 46 had BM (29.7%). A carcinoid syndrome (OR 2.96, p = 0.009) and ≥ 3 extra-skeletal metastatic organs (OR 4.99, p = 0.002) were independently associated with the presence of BM. BM were mainly osteoblastic (78%), rarely symptomatic (8.9%), and had a short-term therapeutic impact for 3 patients (6.5%). The presence of BM (HR 2.67, p = 0.034), older age (HR 1.07, p = 0.016), and higher Ki67 (HR 1.09, p = 0.025) were independent prognostic factors. Unlike poly-BM (HR 1.92, p = 0.007), oligo-BM was not a poor prognosis factor (HR 0.77, p = 0.699) compared to the group without BM. CONCLUSION: 18F-DOPA PET frequently detects BM in patients with metastatic midgut NET. BM have a negative prognostic impact, especially poly-BM. Conversely, oligo-BM do not influence the prognosis and may not impact therapeutic decisions. KEY POINTS: • 18F-DOPA PET detected bone metastases in 46 (29.7%) of 155 patients with metastatic midgut neuroendocrine tumors. • Bone metastases have a negative prognostic impact in metastatic midgut neuroendocrine tumors. • Bone oligo-metastases (< 5) do not influence the prognosis and may not impact therapeutic decisions.

Relationship of Tumor Radiation-absorbed Dose to Survival and Response in Hepatocellular Carcinoma Treated with Transarterial Radioembolization with 90Y in the SARAH Study.

Background Little is known about factors that influence the efficacy of transarterial radioembolization (TARE). Purpose To determine the relationship between tumor radiation-absorbed dose and survival and tumor response in locally advanced inoperable hepatocellular carcinoma treated with TARE. Materials and Methods This was a secondary analysis of prospectively acquired data (between December 2011 and March 2015) from participants who received TARE in the Sorafenib versus Radioembolization in Advanced Hepatocellular Carcinoma (SARAH) trial (ClinicalTrials.gov identifier: NCT01482442). Tumor-absorbed dose was computed using technetium 99m (99mTc) macroaggregated human albumin (MAA) SPECT/CT. Visual agreement among CT, 99mTc-MAA SPECT/CT, and yttrium 90 (90Y) SPECT/CT or PET/CT was scored as optimal, suboptimal, or not optimal. Overall survival (OS) and tumor response at 6-month follow-up CT (Response Evaluation Criteria in Solid Tumors, version 1.1) were assessed. OS was evaluated using Kaplan-Meier tests. A propensity score comparing participants receiving a tumor dose greater than or equal to 100 Gy (best cut-off according to the receiver operating characteristic curve and median tumor radiation-absorbed dose values in the study groups) with those receiving sorafenib was calculated. Results One hundred twenty-one participants (median age, 67 years; interquartile range [IQR]: 61-73 years; 110 men) were evaluated in the dose-survival group, and 109 (median age, 66 years; IQR: 61-71 years; 100 men) were evaluated in the dose-tumor response group. In the dose-survival group, median OS was 9.3 months (95% confidence interval [CI]: 6.7 months, 10.7 months), and median tumor radiation-absorbed dose was 112 Gy (IQR: 68-220 Gy). Participants who received at least 100 Gy (n = 67) had longer survival than those who received less than 100 Gy (median, 14.1 months [95% CI: 9.6 months, 18.6 months] vs 6.1 months [95% CI: 4.9 months, 6.8 months], respectively; P < .001), and those with optimal agreement (n = 24) had the longest median OS (24.9 months; 95% CI: 9.6 months, 33.9 months). In the dose-tumor response group, tumor radiation-absorbed dose was higher in participants with disease control versus those with progressive disease (median, 121 Gy [IQR: 86-190 Gy] vs 85 Gy [IQR: 58-164 Gy]; P = .02). The highest disease control rate was observed in 31 of 40 participants (78%) with a tumor radiation-absorbed dose greater than or equal to 100 Gy and optimal agreement. Conclusion Higher tumor radiation-absorbed dose computed at technetium 99m macroaggregated human albumin SPECT/CT was associated with better overall survival and disease control in hepatocellular carcinoma treated with transarterial radioembolization with yttrium 90 in the Sorafenib versus Radioembolization in Advanced Hepatocellular Carcinoma trial. © RSNA, 2020 Online supplemental material is available for this article. See also the editorial by Sofocleous and Kamarinos in this issue.

[18F]FDG PET/CT predicts progression-free survival in patients with idiopathic pulmonary fibrosis.

BACKGROUND: Idiopathic pulmonary fibrosis (IPF) is a devastating disease characterized by an unpredictable course. Prognostic markers and disease activity markers are needed. The purpose of this single-center retrospective study was to evaluate the prognostic value of lung fluorodeoxyglucose ([18F]-FDG) uptake assessed by standardized uptake value (SUV), metabolic lung volume (MLV) and total lesion glycolysis (TLG) in patients with IPF. METHODS: We included 27 IPF patients (IPF group) and 15 patients with a gastrointestinal neuroendocrine tumor without thoracic involvement (control group). We quantified lung SUV mean and SUV max, MLV and TLG and assessed clinical data, high-resolution CT (HRCT) fibrosis and ground-glass score; lung function; gender, age, physiology (GAP) stage at inclusion and during follow-up; and survival. RESULTS: Lung SUV mean and SUV max were higher in IPF patients than controls (p <0.00001). For patients with IPF, SUV mean, SUV max, MLV and TLG were correlated with severity of lung involvement as measured by a decline in forced vital capacity (FVC) and diffusing capacity of the lungs for carbon monoxide (DLCO) and increased GAP score. In a univariate and in a multivariate Cox proportional-hazards model, risk of death was increased although not significantly with high SUV mean. On univariate analysis, risk of death was significantly associated with high TLG and MLV, which disappeared after adjustment functional variables or GAP index. Increased MLV and TLG were independent predictors of death or disease progression during the 12 months after PET scan completion (for every 100-point increase in TLG, hazard ratio [HR]: 1.11 (95% CI 1.06; 1.36), p = 0.003; for every 100-point increase in MLV, HR: 1.20 (1.04; 1.19), p = 0.002). On multivariable analysis including TLG or MLV with age, FVC, and DLCO or GAP index, TLG and MLV remained associated with progression-free survival (HR: 1.1 [1.03; 1.22], p = 0.01; and 1.13 [1.0; 1.2], p = 0.005). CONCLUSION: FDG lung uptake may be a marker of IPF severity and predict progression-free survival for patients with IPF.

A gate evaluation of the sources of error in quantitative 90Y PET.

PURPOSE: Accurate reconstruction of the dose delivered by 90Y microspheres using a postembolization PET scan would permit the establishment of more accurate dose-response relationships for treatment of hepatocellular carcinoma with 90Y. However, the quality of the PET data obtained is compromised by several factors, including poor count statistics and a very high random fraction. This work uses Monte Carlo simulations to investigate what impact factors other than low count statistics have on the quantification of 90Y PET. METHODS: PET acquisitions of two phantoms-a NEMA PET phantom and the NEMA IEC PET body phantom-containing either 90Y or 18F were simulated using gate. Simulated projections were created with subsets of the simulation data allowing the contributions of random, scatter, and LSO background to be independently evaluated. The simulated projections were reconstructed using the commercial software for the simulated scanner, and the quantitative accuracy of the reconstruction and the contrast recovery of the reconstructed images were evaluated. RESULTS: The quantitative accuracy of the 90Y reconstructions were not strongly influenced by the high random fraction present in the projection data, and the activity concentration was recovered to within 5% of the known value. The contrast recovery measured for simulated 90Y data was slightly poorer than that for simulated 18F data with similar count statistics. However, the degradation was not strongly linked to any particular factor. Using a more restricted energy range to reduce the random fraction in the projections had no significant effect. CONCLUSIONS: Simulations of 90Y PET confirm that quantitative 90Y is achievable with the same approach as that used for 18F, and that there is likely very little margin for improvement by attempting to model aspects unique to 90Y, such as the much higher random fraction or the presence of bremsstrahlung in the singles data.

A gate evaluation of the sources of error in quantitative90 Y PET.

PURPOSE: Accurate reconstruction of the dose delivered by 90 Y microspheres using a postembolization PET scan would permit the establishment of more accurate dose-response relationships for treatment of hepatocellular carcinoma with 90 Y. However, the quality of the PET data obtained is compromised by several factors, including poor count statistics and a very high random fraction. This work uses Monte Carlo simulations to investigate what impact factors other than low count statistics have on the quantification of90 Y PET. METHODS: PET acquisitions of two phantoms-a NEMA PET phantom and the NEMA IEC PET body phantom-containing either 90 Y or 18 F were simulated using gate. Simulated projections were created with subsets of the simulation data allowing the contributions of random, scatter, and LSO background to be independently evaluated. The simulated projections were reconstructed using the commercial software for the simulated scanner, and the quantitative accuracy of the reconstruction and the contrast recovery of the reconstructed images were evaluated. RESULTS: The quantitative accuracy of the 90 Y reconstructions were not strongly influenced by the high random fraction present in the projection data, and the activity concentration was recovered to within 5% of the known value. The contrast recovery measured for simulated 90 Y data was slightly poorer than that for simulated 18 F data with similar count statistics. However, the degradation was not strongly linked to any particular factor. Using a more restricted energy range to reduce the random fraction in the projections had no significant effect. CONCLUSIONS: Simulations of 90 Y PET confirm that quantitative 90 Y is achievable with the same approach as that used for 18 F, and that there is likely very little margin for improvement by attempting to model aspects unique to 90 Y, such as the much higher random fraction or the presence of bremsstrahlung in the singles data.

Implementation and validation of collapsed cone superposition for radiopharmaceutical dosimetry of photon emitters.

Two collapsed cone (CC) superposition algorithms have been implemented for radiopharmaceutical dosimetry of photon emitters. The straight CC (SCC) superposition method uses a water energy deposition kernel (EDKw) for each electron, positron and photon components, while the primary and scatter CC (PSCC) superposition method uses different EDKw for primary and once-scattered photons. PSCC was implemented only for photons originating from the nucleus, precluding its application to positron emitters. EDKw are linearly scaled by radiological distance, taking into account tissue density heterogeneities. The implementation was tested on 100, 300 and 600 keV mono-energetic photons and (18)F, (99m)Tc, (131)I and (177)Lu. The kernels were generated using the Monte Carlo codes MCNP and EGSnrc. The validation was performed on 6 phantoms representing interfaces between soft-tissues, lung and bone. The figures of merit were γ (3%, 3 mm) and γ (5%, 5 mm) criterions corresponding to the computation comparison on 80 absorbed doses (AD) points per phantom between Monte Carlo simulations and CC algorithms. PSCC gave better results than SCC for the lowest photon energy (100 keV). For the 3 isotopes computed with PSCC, the percentage of AD points satisfying the γ (5%, 5 mm) criterion was always over 99%. A still good but worse result was found with SCC, since at least 97% of AD-values verified the γ (5%, 5 mm) criterion, except a value of 57% for the (99m)Tc with the lung/bone interface. The CC superposition method for radiopharmaceutical dosimetry is a good alternative to Monte Carlo simulations while reducing computation complexity.

Medical physics personnel for medical imaging: requirements, conditions of involvement and staffing levels-French recommendations.

The French regulations concerning the involvement of medical physicists in medical imaging procedures are relatively vague. In May 2013, the ASN and the SFPM issued recommendations regarding Medical Physics Personnel for Medical Imaging: Requirements, Conditions of Involvement and Staffing Levels. In these recommendations, the various areas of activity of medical physicists in radiology and nuclear medicine have been identified and described, and the time required to perform each task has been evaluated. Criteria for defining medical physics staffing levels are thus proposed. These criteria are defined according to the technical platform, the procedures and techniques practised on it, the number of patients treated and the number of persons in the medical and paramedical teams requiring periodic training. The result of this work is an aid available to each medical establishment to determine their own needs in terms of medical physics.

A new approach for dose calculation in targeted radionuclide therapy (TRT) based on collapsed cone superposition: validation with (90)Y.

To speed-up the absorbed dose (AD) computation while accounting for tissue heterogeneities, a Collapsed Cone (CC) superposition algorithm was developed and validated for (90)Y. The superposition was implemented with an Energy Deposition Kernel scaled with the radiological distance, along with CC acceleration. The validation relative to Monte Carlo simulations was performed on 6 phantoms involving soft tissue, lung and bone, a radioembolisation treatment and a simulated bone metastasis treatment. As a figure of merit, the relative AD difference (ΔAD) in low gradient regions (LGR), distance to agreement (DTA) in high gradient regions and the γ(1%,1 mm) criterion were used for the phantoms. Mean organ doses and γ(3%,3 mm) were used for the patient data. For the semi-infinite sources, ΔAD in LGR was below 1%. DTA was below 0.6 mm. All profiles verified the γ(1%,1 mm) criterion. For both clinical cases, mean doses differed by less than 1% for the considered organs and all profiles verified the γ(3%,3 mm). The calculation time was below 4 min on a single processor for CC superposition and 40 h on a 40 nodes cluster for MCNP (10(8) histories). Our results show that the CC superposition is a very promising alternative to MC for (90)Y dosimetry, while significantly reducing computation time.

Study of the impact of tissue density heterogeneities on 3-dimensional abdominal dosimetry: comparison between dose kernel convolution and direct Monte Carlo methods.

UNLABELLED: Dose kernel convolution (DK) methods have been proposed to speed up absorbed dose calculations in molecular radionuclide therapy. Our aim was to evaluate the impact of tissue density heterogeneities (TDH) on dosimetry when using a DK method and to propose a simple density-correction method. METHODS: This study has been conducted on 3 clinical cases: case 1, non-Hodgkin lymphoma treated with (131)I-tositumomab; case 2, a neuroendocrine tumor treatment simulated with (177)Lu-peptides; and case 3, hepatocellular carcinoma treated with (90)Y-microspheres. Absorbed dose calculations were performed using a direct Monte Carlo approach accounting for TDH (3D-RD), and a DK approach (VoxelDose, or VD). For each individual voxel, the VD absorbed dose, D(VD), calculated assuming uniform density, was corrected for density, giving D(VDd). The average 3D-RD absorbed dose values, D(3DRD), were compared with D(VD) and D(VDd), using the relative difference Δ(VD/3DRD). At the voxel level, density-binned Δ(VD/3DRD) and Δ(VDd/3DRD) were plotted against ρ and fitted with a linear regression. RESULTS: The D(VD) calculations showed a good agreement with D(3DRD). Δ(VD/3DRD) was less than 3.5%, except for the tumor of case 1 (5.9%) and the renal cortex of case 2 (5.6%). At the voxel level, the Δ(VD/3DRD) range was 0%-14% for cases 1 and 2, and -3% to 7% for case 3. All 3 cases showed a linear relationship between voxel bin-averaged Δ(VD/3DRD) and density, ρ: case 1 (Δ = -0.56ρ + 0.62, R(2) = 0.93), case 2 (Δ = -0.91ρ + 0.96, R(2) = 0.99), and case 3 (Δ = -0.69ρ + 0.72, R(2) = 0.91). The density correction improved the agreement of the DK method with the Monte Carlo approach (Δ(VDd/3DRD) < 1.1%), but with a lesser extent for the tumor of case 1 (3.1%). At the voxel level, the Δ(VDd/3DRD) range decreased for the 3 clinical cases (case 1, -1% to 4%; case 2, -0.5% to 1.5%, and -1.5% to 2%). No more linear regression existed for cases 2 and 3, contrary to case 1 (Δ = 0.41ρ - 0.38, R(2) = 0.88) although the slope in case 1 was less pronounced. CONCLUSION: This study shows a small influence of TDH in the abdominal region for 3 representative clinical cases. A simple density-correction method was proposed and improved the comparison in the absorbed dose calculations when using our voxel S value implementation.

Clinical feasibility of fast 3-dimensional dosimetry of the liver for treatment planning of hepatocellular carcinoma with 90Y-microspheres.

UNLABELLED: Several treatment strategies are used for selective internal radiation therapy with (90)Y-microspheres. The diversity of approaches does not favor the standardization of the prescribed activity calculation. To this aim, a fast 3-dimensional (3D) dosimetry method was developed for (90)Y-microsphere treatment planning and was clinically evaluated retrospectively. METHODS: Our 3D approach is based on voxel S values (VSVs) and has been implemented in the software tool VoxelDose. VSVs were previously calculated at a fine voxel size. The time-integrated activity (TIA) map is derived from pretherapeutic (99m)Tc-macroaggregated-albumin SPECT/CT. The fine VSV map is resampled at the voxel size of the TIA map. Then, the TIA map is convolved with the resampled VSV map to construct the 3D dose map. Data for 10 patients with 12 tumor sites treated by (90)Y-microspheres for hepatocellular carcinoma were collected retrospectively. 3D dose maps were computed for each patient, and tumoral liver and nontumoral liver (TL and NTL, respectively) were delineated, allowing the computation of descriptive statistics (i.e., mean absorbed dose, minimum absorbed dose, and maximum absorbed dose) and dose-volume histograms. Mean absorbed doses in TL and NTL from VoxelDose were compared with those calculated with the standard partition model. RESULTS: The estimated processing time for a complete 3D dosimetry calculation is on the order of 15 min, including 10 s for the dose calculation (i.e., VSV resampling and convolution). An additional 45 min was needed for the semiautomatic and manual segmentation of TL and NTL. The mean absorbed dose (±SD) was 422 ± 263 Gy for TL and 50.1 ± 36.0 Gy for NTL. The comparison between VoxelDose and partition model shows a mean relative difference of 1.5% for TL and 4.4% for NTL. Results show a wide spread of voxel-dose values around mean absorbed dose. The minimum absorbed dose within TL ranges from 32 to 267 Gy (n = 12). The fraction of NTL volume irradiated with at least 80 Gy ranges from 4% to 70% (n = 10), and the absorbed dose from which 25% of NTL was the least irradiated ranges from 14 to 178 Gy. CONCLUSION: This article demonstrates the feasibility of a fast 3D dosimetry method for (90)Y-microspheres and highlights the potential value of a 3D treatment planning strategy.

Fine-resolution voxel S values for constructing absorbed dose distributions at variable voxel size.

UNLABELLED: This article presents a revised voxel S values (VSVs) approach for dosimetry in targeted radiotherapy, allowing dose calculation for any voxel size and shape of a given SPECT or PET dataset. This approach represents an update to the methodology presented in MIRD pamphlet no. 17. METHODS: VSVs were generated in soft tissue with a fine spatial sampling using the Monte Carlo (MC) code MCNPX for particle emissions of 9 radionuclides: (18)F, (90)Y, (99m)Tc, (111)In, (123)I, (131)I, (177)Lu, (186)Re, and (201)Tl. A specific resampling algorithm was developed to compute VSVs for desired voxel dimensions. The dose calculation was performed by convolution via a fast Hartley transform. The fine VSVs were calculated for cubic voxels of 0.5 mm for electrons and 1.0 mm for photons. Validation studies were done for (90)Y and (131)I VSV sets by comparing the revised VSV approach to direct MC simulations. The first comparison included 20 spheres with different voxel sizes (3.8-7.7 mm) and radii (4-64 voxels) and the second comparison a hepatic tumor with cubic voxels of 3.8 mm. MC simulations were done with MCNPX for both. The third comparison was performed on 2 clinical patients with the 3D-RD (3-Dimensional Radiobiologic Dosimetry) software using the EGSnrc (Electron Gamma Shower National Research Council Canada)-based MC implementation, assuming a homogeneous tissue-density distribution. RESULTS: For the sphere model study, the mean relative difference in the average absorbed dose was 0.20% ± 0.41% for (90)Y and -0.36% ± 0.51% for (131)I (n = 20). For the hepatic tumor, the difference in the average absorbed dose to tumor was 0.33% for (90)Y and -0.61% for (131)I and the difference in average absorbed dose to the liver was 0.25% for (90)Y and -1.35% for (131)I. The comparison with the 3D-RD software showed an average voxel-to-voxel dose ratio between 0.991 and 0.996. The calculation time was below 10 s with the VSV approach and 50 and 15 h with 3D-RD for the 2 clinical patients. CONCLUSION: This new VSV approach enables the calculation of absorbed dose based on a SPECT or PET cumulated activity map, with good agreement with direct MC methods, in a faster and more clinically compatible manner.

Spatial gradients of blood vessels and hematopoietic stem and progenitor cells within the marrow cavities of the human skeleton.

This report evaluates the spatial profile of blood vessel fragments (BVFs) and CD34(+) and CD117(+) hematopoietic stem and progenitor cells (HSPCs) in human cancellous bone. Bone specimens were sectioned, immunostained (anti-CD34 and anti-CD117), and digitally imaged. Immunoreactive cells and vessels were then optically and morphometrically identified and labeled on the corresponding digital image. The distance of each BVF, or CD34(+) or CD117(+) HSPC to the nearest trabecular surface was measured and binned in 50-microm increments. The relative concentration of HSPCs and BVFs within cancellous marrow was observed to diminish with increasing distance in the marrow space. On average, 50% of the CD34(+) HSPC population, 60% of the CD117(+) HSPC population, and 72% of the BVFs were found within 100 microm of the bone surfaces. HSPCs were also found to exist in close proximity to BVFs, which supports the notion of a shared HSPC and vessel spatial niche.

Comparison between 2D and 3D dosimetry protocols in 90Y-ibritumomab tiuxetan radioimmunotherapy of patients with non-Hodgkin's lymphoma.

UNLABELLED: We compared the radiation-absorbed dose obtained from a two dimensional (2D) protocol, based on planar whole-body (WB) scans and fixed reference organ masses with dose estimates, using a 3D single-photon emission computed tomography (SPECT) imaging protocol and patient-specific organ masses. METHODS: Six (6) patients with follicular non-Hodgkin's lymphoma underwent a computed tomography (CT) scan, 5 2D planar WB, and 5 SPECT scans between days 0 and 6 after the injection of (111)In-ibritumomab tiuxetan. The activity values in the liver, spleen, and kidneys were calculated from the 2D WB scans, and also from the 3D SPECT images reconstructed, using quantitative image processing. Absorbed doses after the administration of (90)Y-ibritumomab tiuxetan were calculated from the (111)In WB activity values combined with reference organ masses and also from the SPECT activity values and organ masses as estimated from the patient CT scan. To assess the quantitative accuracy of the WB and SPECT scans, an abdominal phantom was also studied. RESULTS: The differences between organ masses estimated from the patient CT and from the reference MIRD models were between -10% and +98%. Using the phantom, errors in organ and tumor activity estimates were between -86% and 10% for the WB protocol and between -43% and -6% for the SPECT protocol. Patient liver, spleen, and kidney activity values determined from SPECT were systematically less than those from the WB scans. Radiation-absorbed doses calculated with the 3D protocol were systematically lower than those calculated from the WB protocol (29%+/-26%, 73%+/-26%, and 33%+/-53% differences in the liver, spleen, and kidney, respectively), except in the kidneys of 2 patients and in the liver of 1 patient. CONCLUSIONS: Accounting for patient-specific organ mass and using SPECT activity quantification have both a great impact on estimated absorbed doses.

Semiautomated thoracic and abdominal computed tomography segmentation using the belief functions theory: application to 3D internal dosimetry.

AIM: Segmentation of computed tomography (CT) images is an important step in three-dimensional (3D) internal dosimetry. To this end, a semiautomated method was developed to delineate organs, using the belief functions theory. MATERIALS AND METHODS: The membership degree of each voxel to each volume of interest is estimated by computing a basic belief assignment (bba). For each voxel V, bbas corresponding to each neighbor are aggregated to obtain a unique bba, using a merging procedure. Before aggregating information, a 3D filter is applied, in order to take into account the fact that the more the voxel V(i) is close to V, the more the information coming from V(i) is reliable. The aim is to weaken the contribution of voxels according to their distance with respect to the voxel to be classified. The algorithm was applied on 10 CT scans (pixel size, 0.98 x 0.98 mm2, slice thickness 3 mm, 120 kV). For each organ (i.e., the lung, liver, kidney, and spleen), the algorithm was applied on a part of the CT volume. First, the lung was segmented using two classes with characteristic values, C, defined by the K-means clustering algorithm. Second, the liver and kidneys were segmented using three classes with C-values defined by local mean Hounsfield Unit (HU) measurements, corresponding to fat, liver, and kidney. Third, the spleen was segmented using three classes corresponding to fat, kidney, and spleen, using local mean HU measurements. The semiautomated segmentation was compared with manual segmentation using the volume difference and the agreement (overlap) index. RESULTS: For organ segmentation, the computation duration was between 5 and 20 minutes (2.5 GHz, RAM of 1 GByte), depending on the number of classes and the volume size to classify. On the 10 patients, manual correction was needed for none on the lung, 1 on the spleen, 2 on the kidneys, and 7 on the liver, mostly owing to intercostal structures. The mean relative volume difference (+/-1 standard deviation [SD]) between manual and automated segmentation was 5.0% +/- 3.7%, 5.7% +/- 3.6%, 6.2% +/- 3.3%, and 7.5% +/- 6.5% for the lung, liver, kidney, and spleen, respectively. The corresponding mean agreement index (+/-1 SD) was 0.94 +/- 0.03, 0.90 +/- 0.01, 0.88 +/- 0.03, and 0.84 +/- 0.04. CONCLUSIONS: The algorithm allows for the delineating of organs on thoracic and abdominal CT images, and will be integrated in a 3D internal dosimetry dedicated software.